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Abstract. We investigate the properties of nonlinear excitations in different types of soliton bearing systems 
with long-range dispersive interaction. We show that length-scale competition in such systems universally 
results in a multi-component structure of nonlinear excitations and can lead to a new type of multistability: 
coexistence of different nonlinear excitations at the same value of the spectral parameter (i.e., velocity in 
the case of anharmonic lattices or frequency in nonlinear Schrodinger models). 
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1 Introduction 

Physical systems with competition between discreteness, 
nonlinearity and dispersion are abundant in Nature. Ex- 
amples are nonlinear charge and energy transport in bi- 
ological macromolecules, elastic energy transfer in anhar- 
monic chains, charge transport in hydrogen-bonded sys- 
tems, coupled optical fibers, nonlinear photonic crystals, 
arrays of coupled Josephson junctions. Determination of 
their dynamical properties is of importance because of 
their wide applicability in various physical problems. In 
nonlinear systems with short-ranged (nearest-neighbor) 
dispersive interactions there are basically two types of 
nonlinear excitations. In the large dispersion (weak non- 
linearity) limit when the discreteness of the lattice plays 
no essential role, the balance between nonlinearity and 
dispersion provides the existence of low-energy soliton- 
like excitations. The width of these excitations is large 
comparing with the lattice spacing, they freely propagate 
along the system without energy loss, their collisions are 
almost elastic. When the nonlinearity is strong (dispersion 
is weak) new nonlinear excitations, intrinsically localized 
modes, also called discrete breathers, may appear (see 6.17. 
Ref. [y for a review). Their width is comparable with the 
lattice spacing, they are long-lived and robust. 

The phenomenon of multistability, implying the coex- 
istence of several stationary localized states haviiig the 
same energy, appears m different contexts 
|9|, Pt1[i^ ,1i3|, . For the one-dimensional discrete non- 
linear Schrodinger (NLS) equation with nearest-neighbor 
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interactions (NNI) and arbitrary power of nonlinearity, 
multistability occurs when the power exceeds some criti- 
cal value §11. 

Another mechanism of such multistability comes out 
for the systems with long-range dispersive interactions, 
where a new characteristic length scale appears: the radius 
of the dispersive interaction. As it was shown in Rcfs. §,|^, 
the appearance of the new characteristic length 
scale brings into coexistence two stable stationary excita- 
tions one of which is a broad, continuum-like soliton and 
another is a narrow discrete excitation. This is a generic 
property of nonlinear systems with nonlocal dispersive in- 
teractions. The bistability occurs in such different systems 
as anharmonic chains with nonlocal interatomic interac- 
tion NLS models with long-range coupHng 
and nonlinear waveguides in photonic crystals ||l0|] . In all 
these systems the bistability means that there is an in- 
terval of the energy H where two stable nonlinear ex- 
citations (one of which is broad and another is intrinsi- 
cally localized) coexist for each value of H . The greater 
the radius of the long-range interactions is, the more pro- 
nounced is the bistability phenomenon. A general mecha- 
nism for obtaining a controlled switching between bistable 
localized excitations was proposed in Refs. 1 1 1 , 1^ . Impor- 
tantly, very similar types of excitations may appear in pe- 
riodic nonlinear dielectric superlattices (two-dimensional 
nonlinear Kronig-Pcnney model p3[|) where the nonlocal- 
ity arises as a result of a reduction of the original two- 
dimensional problem to an effective one-dimensional de- 
scription. Bistability in the spectrum of nonlinear local- 
ized excitations in magnets with a weak exchange inter- 
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action and strong magnetic anisotropy was found in Ref. 

Quite recently coexistence of different nonlinear exci- 
tations at the same value of the soliton velocity has been 
shown to take place in anharmonic chains with ultra-long- 
range interatomic interactions . Specifically, it was shown 
that three different types of pulse solitons can coexist for 
each value of the pulse velocity in a certain velocity in- 
terval: one type is unstable but the two others can be 
stable. It should be emphasized the following intriguing 
aspect of the problem: in contrast to the previously dis- 
cussed cases, now the bistability takes place at the same 
value of a spectral parameter (the soliton velocity in this 
system is mathematically the spectral parameter of a non- 
linear eigenvalue problem). Therefore, the system being 
discussed provides us with the example of a new type of 
multistability whose physical consequences can be quite 
different from the multistability discussed above. Indeed, 
if such a new type of multistability can exist, for instance, 
for nonlinear waveguides in photonic crystals pO| , then 
two different nonlinear modes with the same frequency 
of the light will coexist in a waveguide, which is practi- 
cally more important than to have coexistence of nonlinear 
modes with the same energies but different frequencies. 

The aim of the present paper is to demonstrate that 
the new type of multistability is a fundamental property 
of nonlinear systems with ultra-long-range dispersion. To 
this end we consider two additional (in comparison with 
Ref. Il^) models: a discrete NLS equation with nonlo- 
cal dispersion and an anharmonic lattice with screened 
Coulomb interaction. We show that in both cases an ultra- 
long-range dispersive interaction results in a multi-valued 
dependence of the energy of stationary states on a spectral 
parameter (frequency in the case of the NLS equation and 
velocity in the case of the anharmonic lattice). 

The outline of the paper is the following. In Sec. 2 we 
present the results of numerical simulations of stationary 
states in the discrete NLS model with a long-range disper- 
sive interaction and in the anharmonic chain with a long- 
range dispersive interaction and cubic nearest-neighbor in- 
teraction. We discuss the appearance of a multi- valued de- 
pendence of the energy of the systems. In Sec. 3 we develop 
an analytical approach to the problem. We use a quasicon- 
tinuum approach which preserves both length scales that 
enter our problem. We show that the low-energy excitation 
can be considered as a compound state: a short-ranged 
breather- like component and a long-ranged component. 
Section 4 presents the concluding discussion. 



2 Models and numerical results 

We consider in parallel two nonlinear models which can 
bear pulse- like excitations: the Boussinesq equation model 
which describes the properties of anharmonic molecular 
chains and the nonlinear Schrodinger equation model which 
represents generic properties of a system of nonlinear os- 
cillators. 



2.1 Anharmonic chains with ultra-long-range 
dispersive interactions 

Let us consider a chain of equally spaced particles of unit 
mass whose displacements from equilibrium are m„ (t) and 
the equilibrium spacings are unity. The Hamiltonian of 
the system is given by 

n ^ 

1 ^ 

+ ^ / ^ Jm-n( 



with the anharmonic interaction 

2 /o „,.3 , 



V{w) = w72- w73 



(1) 



(2) 



between nearest neighbors and the harmonic long-range 
interaction (LRI) 



Jm-n{a, s) = J(a, s) 



^-a\m-n\ 

\m — nY 



between all particles of the chain. Here 



J(a,s) = J/F{e-°',s) 



where 



F{z,s) 



(3) 



(4) 



(5) 



is the so-called Joncqiere's function which properties are 
described in Ref. . J characterizes the intensity of the 
LRI's whereas a and s determine their inverse radius. The 
coefficient J(a, s) is chosen in the form (^ to have the sum 
X^ra Jnioi,s) independent of a and s. The parameters a 
and s are introduced to cover different physical situations 
from the limit of nearest-neighbor interactions (a 3> 1 or 
s ^ 1) to the limit of ultra-long-range interactions (a <C 1 
and s < 3). The Hamiltonian (P generates equations of 
motion of the form 



d^ 



^JL + 2WiWn) - W{Wn+l) - Wiw^.i) 
+ ^ Jm,n{Wn ~ Wm) = , 



(6) 



are relative displacements and 

2 



where Wn = Un+l — Ur, 

W{w) = dV{w)/dw = w — w'^ . We use the quasicon- 
tinuum approach, regarding n as a continuous variable: 
X, Wn{t) w{x,t), and using the operator identities 



w{x -\- m,t) = e 
Q{a,s,dx) = 



w{x, t) 



m— 1 



1 — cosh mc^T 



= 1 - 



F{e- 



2F(e-",s) 



(7) 
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Fig. 1. Energy of the pulse solitons versus velocity in the 
anharmonic lattice with Kac-Baker (s — 0) long-range inter- 
action found numerically for different values of a and J (the 
value of J was chosen to get constant c = 1.515): a = 0.05 
and J = 0.0016 (full line); a = 0.1 and J = 0.0062 (dashed 
line); a = 0.17 and J — 0.0172 (long-dashed line); a = 0.3 and 
J — 0.05 (dot-dashed line). 



1.5r 




00 0.05 0.10 0.15 0.20 0.25 



2 2 
V -C 

Fig. 2. Energy of the pulse solitons versus velocity in the 
anharmonic lattice with screened Coulomb (s = 3) interaction 
for J = 0.05 and different values of a: a = 0.001 (full line), 
a — 0.005 (dashed line), and a — oo (long-dashed line). 



In this approach the equation of motion (|g) can be written 
in the pseudo-differential form 

[df - JQ{a,s,d:,)]w{x,t) - ismh" (^^^W{w) ^ , (8) 

where dx and dt are the derivatives with respect to x and 
respectively. Eq. (^ can be called non-local Boussinesq 
equation. 

The energy (|l|) together with the momentum 



P = 



(9) 



are conserved quantities. We are interested in the station- 
ary soliton solutions w{x, t) = w{x — vt) propagating with 



Fig. 3. Shapes of the low-energy (full line), unstable 
intermediate-energy (dotted line), and high-energy (dashed 
line) pulse solitons which coexist at the same velocity v — 1.64 
in the anharmonic lattice with Kac-Baker (s — 0) long-range 
interaction for a — 0.05 and J = 0.0016 (see Fig. hj). 
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Fig. 4. Shapes of the low-energy (full line) and high-energy 
(dashed line) pulse solitons which coexist at the same veloc- 
ity v'^ — = 0.031 in the anharmonic lattice with screened 
Coulomb (s = 3) interaction for J = 0.05 and a = 0.001 (see 
Fig- I)- 



velocity v. Substituting z = x 
in the form 



vt we can write Eq. (0) 



[v'^dl ~ JQ{a, s,d;)]w{z)-A sinh^ (^^^W{w) = Q . (10) 

In this way we reduce our problem to a nonlinear eigen- 
value problem with v being a spectral parameter. 

We consider two physically important cases: screened 
Coulomb interactions (s = 3) and the Kac-Baker LRI's 
(s = 0). The speed of sound c (which is an upper limit 
of the group velocity of linear waves), determined by the 
expression = I + JF (e~", s — 2) /F {e~", s), grows in- 
definitely in both cases as a decreases. Using the Green's 
function method Q one can show that spatially localized 
stationary solutions (solitons) exist only for supersonic ve- 
locities V > c. 

We have numerically integrated Eq. ( [l0| ) using the 
method developed in Ref . . Figures [l] and || show the 
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dependence of the energy H{v) (the dependence Pjv) has 
the same form) of the soliton solutions of Eq. (Jl^) ob- 
tained for the Kac-Baker interaction (s = 0) and screened 
Coulomb interaction (s = 3), respectively. It is seen that 
the behavior in both cases is qualitatively the same. The 
soliton energy grows monotonically with the velocity in the 
case of large a (see, e.g., a = 0.3 in Fig. |l|). In this case the 
soliton properties are qualitatively the same as in the limit 
of nearest-neighbor interactions (NNI's) for which Eq. (||) 
reduces to the Boussinesq equation. It is well known that 
this equation has a sech-shaped soliton solution w{z) — 

— 1.5(u^ — c^)/ cosh^((Tz), where a — \/'i{v'^ — c^) is the 
inverse soliton width. As indicated above, the energy H ^ 
{v^ — c^)'^/^ of these solitons is a monotonic function of 
the velocity, which means that there is only one soliton 
solution for each given value of energy or velocity. 

The situation changes when a < Ucr (e.g., for the 
Kac-Baker interaction the critical value of inverse radius 
Oicr{J) is determined by Eq. (9) in Ref. jl^). In this case, 
as was shown first in Ref. the dependence H{v) be- 
comes non-monotonic (of A^-shape, with a local maximum 
and a local minimum) , so that there exist two branches of 
stable supersonic solitons: low- velocity and high- velocity 
solitons. These two soliton branches are separated by a 
gap (an interval of velocities where H(v) decreases) with 
unstable soliton states. It is interesting that because of 
the non-monotonic dependence of H(v) there exist low- 
velocity and high- velocity solitons of equal energy. Further 
increasing of the radius of the dispersive interaction (de- 
creasing a) leads to a Z-shape dependence of H{v) (see 
Figures 0and|). In this case the gap which separated two 
above-mentioned soliton branches collapses; instead, there 
appears an interval of velocities where three soliton solu- 
tions with different energies and shapes (see Figures ^ and 
^) coexist for each value of velocity. This phenomenon was 
first observed in Ref. ||l5| for the Kac-Baker dispersive in- 
teraction. Here we show that the same phenomenon takes 
place also in the case of the screened Coulomb interaction, 
but the radius of the interaction should be much larger: 
e.g., for the Coulomb dispersive interaction with the in- 
tensity J = 0.05 the gap vanishes for a < 0.0015 while 
for the Kac-Baker interaction with the same intensity J it 
happens for a < 0.15. 

2.2 Nonlinear Schrodinger model with ultra-long-range 
dispersive interaction 

Let us now consider a system of coupled nonlinear oscil- 
lators which is described by the Hamiltonian 



H = 




2 ^ 



,(a,s)|V„-^™p-|V'nn ,(11) 



where n and m are site indices, and ■(/;„ is the excitation 
wave function. The first term in Eq. (^ij) describes a dis- 
persive interaction between nonlinear oscillators. The sec- 
ond one represents a nonlinear interaction in the system. 
The case a = in the dispersive interaction has been stud- 
ied in detail in Ref. M. Here we consider the dispersive 
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Fig. 5. Energy H of the stationary solitons versus nonlinear 
frequency A in the NLS model with Kac-Baker (s = 0) long- 
range interaction for different values of a: 0.1 (full line), 0.3 
(dashed line), 0.7 (long-dashed line), and oo (dot-dashed line). 
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Fig. 6. Number of quanta N of the stationary solitons versus 
nonlinear frequency A in the NLS model with Kac-Baker (s = 
0) long-range interaction for different values of a: 0.1 (full line), 
0.3 (dashed line), 0.7 (long-dashed line), and oo (dot-dashed 
line) . 



interaction J„_„i(a, s) in the Kac-Baker form (s = 0): 

J„_„ = J(e"-l)e-"l"^"l . (12) 

It is important to investigate this type of the dispersive 
interaction because real physical systems which are de- 
scribed by the Hamiltonian (11) with the dispersive in- 
teraction given by Eq. ([l^ ) are not rare. For example, 
such type of excitation transfer may take place in systems 
where the dispersion curves of two elementary excitations 
are close or intersect and effective long-range transfer oc- 
curs at the cost of the coupling between excitations. This 
is the case for excitons and photons in semiconductors 
and molecular crystals. Dynamics of nonlinear photonic 
band-gap materials and periodic nonlinear dielectric su- 
perlattices is described [|3| by equations which are closely 
related to the equations which govern the dynamics of the 
discrete NLS model with the dispersive interaction in the 
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Kac-Baker form (p2|). It was shown in Ref. jl^ that the ra- 
dius of the effective dispersive interaction is proportional 
to the distance between nonlinear layers. Therefore one 
can easily tune the range of the dispersive interaction by 
preparing nonlinear lattices with different spacing between 
nonlinear layers. 

In the Hamiltonian (|llf) ipn and ii/jn are canonically 



conjugated variables and the equation of motion i^V'; 
dH/dip^ has the form 



dt 



^n) + |^„|Vn = 0. (13) 



The energy of the system H together with the number of 
quanta 



N 



(14) 



are conserved quantities. 

We study stationary states of the system 



iAt 



(15) 



where A is the nonlinear frequency and (j)n{A) is the am- 
plitude at the n-th site. For small and intermediate ra- 
dius {oi~^) of the dispersive interaction the problem was 
studied in Refs. |ll]. It was shown that for a < 1.7 the 
NLS model exhibits bistability in the spectrum of non- 
linear stationary states: for each value of N in the inter- 
val [-/Vi(a), Nu{a)] there exist three stationary states with 
frequencies yli(iV) < A2{N) < A3{N). Since for discrete 
NLS equations a necessary and sufficient condition for lin- 
ear stability of stationary states with a single maximum 
at a lattice site is dN/dA > (see, e.g., Ref. [H), the 
low-frequency Ai{N) and high-frequency A3{N) states, 
for which dN/dA > 0, are stable (in what follows we 
shall refer to them as to the 'low-frequency branch' and 
'high-frequency branch'). There is an interval ('gap') of 
the nonlinear frequency A which separates two branches 
of nonlinear excitations. Figures ^ and ^ show H versus A 
and N versus A, respectively, obtained from direct numer- 
ical solution of Eq. (^3|) for different values of the inverse 
radius of dispersive interaction a. It is seen that similar 
to the case of the nonlocal Boussinesq equation the gap 
decreases when the radius of the dispersive interaction in- 
creases and eventually it vanishes. 

Comparing Figs. |^ and ^ obtained for the Boussinesq 
model with Fig. || obtained for the NLS model we see 
that several new but common for both models features 
arise as a consequence of the ultra-long-range dispersive 
interaction: 

— The gap which separates two branches of stable sta- 
tionary solutions decreases when the radius of the dis- 
persive interaction increases. 

— For small enough a (in the NLS model for a < 0.57) 
three different stationary states coexist for each value 
of spectral parameter (velocity v in the case of the 
Boussinesq equation and nonlinear frequency A in the 
NLS model). These three states are characterized by 
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Fig. 7. Shapes of the high-energy (full line) and low-energy 
(dashed line) stationary solitons which coexist at the same non- 
linear frequency A = 0.25 in the NLS model with Kac-Baker 
long-range interaction for a — 0.3 (see Fig. H). 



different values of the integrals of motion: the energy 
H (in both models) and either the number of quanta 
N (for the NLS equation) or the momentum P (for the 
Boussinesq equation) . 

The shapes of these states differ significantly (see Figs. 
||-Qand ^). The high-energy state is broad and continu- 
um-like while the low-energy states consist of two com- 
ponents, short-range and long-range ones. The width 
of the short-range component is of a few lattice spacing 
while the long-range component has tails which spread 
out a few tens of lattice constants. 
The state with an intermediate energy is always un- 
stable while the stability of two other states depends 
on the radius of the dispersive interaction a^^ . For ex- 
ample, in the case of the NLS equation for a > 0.03 
the state which belongs to the high-frequency branch 
is unstable {dN/dA < for these a). But this state 
becomes stable when a < 0.03. 



3 Analytical approach 

In our analytical approach we will be mostly concerned 
with the case of the NLS equation model because an ana- 
lytical theory of multi-component solutions of the nonlocal 
Boussinesq equation is partially presented in Refs. lq,0l. 



3.1 Nonlinear Schrodinger equation model 

In the quasicontinuum approach, regarding n as a contin- 
uous variable: n — > a;, ^„(t) 'ip{x,t), Eq. ( p^ can be 
written as 

+ ^ Q{a, s = 0, d^)^P + I VI V = , (16) 



where 



g(a,s = 0,9,) = (e" + l) 



4sinh^(9:r/2) 
k2 -4sinh2(9^/2) 



(17) 
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with K = 2sinh(a/2) being a range parameter, is a hnear 
pseudo-difFercntial operator for the Kac-Baker dispersive 
interaction (p^). 

In order to obtain an analytic solution of the pseudo- 
differential equation ( p^ we use an approach akin to the 
previously proposed in Ref. we multiply Eq. by 

the operator ^ 
operator ( 



, , and expand in a Taylor series this 

, ~ 1 ~ T^dl) as well as the operator 

sinh^(9a:/2) in the denominator of the r.h.s. of Eq. (|l7|). 
We stress that our modification of the standard contin- 
uum approximation is crucial in the preserving the second 
length scale that enters our problem. As a result we get a 
pseudo-differential equation 

*(l - + ^(e" + + I^IV = (18) 

where the small term j^d?;.\ip\'^^ was omitted. 
Looking for stationary solutions 



J 



(19) 



where r = -^(e 



l)i, we can write Eq. (18) in the form 

[dl - sl){dl - s-_)cp + ^(a^ - .^)<?3 = , (20) 
with the parameters s± given by 




In the limit of small a 



12- 



(21) 



(22) 



- A+l' + A 

For small nonlinear frequency A we expect excitations 
with a width much larger than so that only one scale 
sZ.^ should be important. Since the characteristic width 
of the excitations of interest here is large as compared 
we can neglect the operator in the operator 
and get 

1 



with s 



A + 1 



(23) 



This equation was investigated in Ref. It was shown 
that Eq. (p3) has a soliton solution when the inequality 
A = [k^ /s^) > 9 holds. As it is seen from Eq. (p3) in 
the case of small a this inequality reduces to yl < 1/8. 
The number of quanta which corresponds to the solution 
of Eq. (23) has the form 



N = N^ ^ siA 



2V^ 



1 



A~<d 
24 



/(x^^+iKVA + s)^ 

^^(^/I-l)(^/A-3)^ 



(24) 



For large A^ when both length-scales are important, we 
will seek a solution of Eq. (Efl) as the sum 



{(l>s{x) + c^l{x)} 



(25) 



where 4>s(x) is the short-range component and 4>l(x) is 
the long-range one. Here (j^six) will dominate in the cen- 
ter, while (t>L{x) in the tails. 

Inserting Eq. (Esh into Eq. (EG) yields 



(dl - si) {{dl - 4)(0s + 0l) + (05 + M'} 

~iK,^-sl){cbs + cbLf ^0 . (26) 

Assuming that the function (t>six) satisfies the equation 

(dl - 4)0s + 30l(O)4 + 30i(a;)</)5 + 4 = » > (27) 



we obtain from Eq. ( p6[ ) an equation for (f>L (x) in the form 

^ B {k^ ~ si) S{x) , (28) 

where 



B = 



dx {4(x) + 3,pl{x)<Ps{x) + 3Mx)4'lix)} 



dx(j)s{x) . 



(29) 



To obtain Eq. ( p8| ) we took into account the big difference 
in the short-range scale s^^ and the long-range scale sl^ . 
It permitted us to replace the product (j)L{z)4>s{z) by 
(I>l{^)4>s{z) in Eq. (|27|), to neglect the term d^ipL as com- 



pared with s%_(f>L in the l.h.s. of Eq. 



and to replace the 



function (k^ - s^){0|(a;) + 'i4>'i{x)(j3six) + 3(?!)l(x)(/)|(x)} 
b^ the r.h.s. of Eq. (|2q). Also, Eq. (|27| ) was used in Eq. 

We obtain for the short-range component (psix) the 
expression 



(t^six) 



V2s+{1 - 3a2 



Vl — cosh(\/l — Sa'^s+x) + \f2a 



(30) 



where a = 0l(O)/s+ is the amplitude of the long-range 
component. 

Eq. (^8|) on the half x-axis subject to the vanishing at 
infinity boundary conditions can be presented as 

(4 - 3 0i)5,0L 



4,2__,^if_t^^2 +^2,^4 (31) 

with the boundary condition at the origin in the form 

-2 (4-3 0i(o)) a,(/)i(+o) = {te-s'_)B . (32) 



Substituting Eqs. (^-(|3l|) into Eq. ( |32| ) we obtain that 
the amplitude a of the long-range component is deter- 
mined by the equation 



(A - 1) arccos ^J2a'^l{\ - d^) ' 



(33) 
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where R{a) = y/Aa'^- {A + 3)a'^/2 + l . 

The number of quanta which corresponds to the solu- 
tion (25) is determined by the expression 



N 




When K > Kb, where Kb ~ 0.63 is the maximum value of 
the fimction in the r.h.s. of Eq. (|3|) for A = 9 (yl = 1/8), 
Eq. (^) has only one solution in the interval A < 9 {A > 
1/8). Inserting Eq. ( |3^ ) into Eq. (Q) and combining Eqs. 
( P^ ) and (HJ) we obtain that the number of quanta as a 
function of A has an A^-type shape (see Fig. ^ with a finite 
gap which separates two branches of stable solitons. When 
K < Kb the gap collapses and Eq. (|3|) has two solutions 
for A > 9 {A < 1/8). Taking into account that for A > 
9 {A < 1/8) there exists also a low-frequency solution 
of Eq. (p3[), one can conclude that for k < Kb there is 
an interval of A where three stationary states correspond 
to each value of the nonlinear frequency. In this case the 
curve N{A) has a ^-type shape (see Fig. ||). It is worth 
noticing that the bifurcation value for the inverse radius 
of the dispersive interaction ab ~ 0.61 (2 sinh(Q;b/2) = Kb) 
obtained analytically is in a surprisingly good agreement 
with the value obtained from numerical simulations {ab — 
0.57 ±0.01). 

It is interesting also to note that when a ^ the am- 
plitude of the long-range component tends to zero while 
the short-range component in the low frequency limit be- 
comes (5-function like. This result is also in a good quali- 
tative agreement with the results of numerical simulations 
(see Fig. I). 



The parameter cr+ is finite at all velocities v > c and 
tends to y/l2 for v ^ oo. The parameter (T_ vanishes at 
V = c and tends to k for u — > cx). One can find that in 
the travel-wave approach the Hamiltonian is determined 
by the expressions 



H ^vP- M{v) 



where P is the canonical momentum and 



M(v) =^ - I w^(x)dx 
6 ' 



(38) 



(39) 



oo 



is the nonlinear part of the Hamiltonian. 

Near the speed of sound excitations have a width much 
larger than cr^"^ and can be described by Eq. ( ^ in the 
approximation {d^—a^)w « —a'^w. In this approximation 
the soliton solutions exist in a finite interval of velocities, 
c < V < Vcr — \J (4c^ — l)/3 . The explicit form of the low- 
velocity soliton in given in Ref. Using this expression 
it is straightforward to obtain the energy of the soliton 
belonging to the low-velocity branch in the form 



M<{v) 



H< = vP< - M<{v) 



93312^13 

{A-^)^{A + b) 



A^ -2A + 40 
I{w^) 



(40) 



with 

I{w) = 



-dx 



3.2 Non-local Boussinesq equation model 

Here we also restrict ourselves to the case when the dis- 
persive interaction is of Kac-Baker type. In this case the 
quasicontinuum limit of the nonlocal Boussinesq equation 
( p^ ) can be presented as follows |^ 

{dl - ol){dl - a'_)wiz) = ^{d^^ - K^)w'{z) , (35) 

where the parameters a± play the same role as the pa- 
rameters s± in the previous subsection. They are given 

by 




(36) 



where the speed of sound c in the lattice with the Kac- 
Baker interaction is determined by the expression 



1 + J 



1 



(e-" - 1) 



1/2 



(37) 



A 



In 



^yAn{w) +Aw-{A + 2)/3 



A-iA + 2)/S 



n{x) = \jAx'^ - -{A + 2)x + l , 



(41) 



where A = (k^/ct^) > 4 in the velocity interval c < v < 
Vcr, and 



Wi 



±(^A + 2-^{A-l)(A-4)) 



(42) 



is the amplitude of the low-velocity soliton. 

Far from the speed of sound we will use the same ap- 
proach as in the previous subsection (see also Ref. j^). We 
seek a solution of Eq. (10) as a sum 

w{z) =wsiz) + WLiz) , (43) 

where wsiz) is the short-range component of the strain 
and wl{z) is the long-range one. 

Assuming that the function ws{z) satisfies the equa- 
tion 

/ 24 \ 12 

dlws-[<jl + ^WL{0)\ws-^wl = 0, (44) 
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we obtain from Eq. (pq) an equation for wl{z) in the form 



12 



where 



B = 



dzws{z) 



(45) 



(46) 



To obtain Eqs. (44) and ( |45| ) we took into account the big 
difference in the short-range scale a^^ and the long-range 
one aZ^ and proceeded in the same way as we did in the 
previous subsection. From Eq. (|^ we get the short-range 
component ws in the form 

z«5(z) = -i^>^(l-27)sech2(yT^^) , (47) 

where 7 — 12\wl{0)\ / {va+)^ determines the amplitude 
of the long-range component. It can be found from the 
equation [p| 



3{A- 1) 



7+ x/T^ 



(48) 



Note that there is a misprint in the corresponding equation 
of Ref. 1^: one should omit wl{0) (which plays the same 
role as 7 in this paper) in the nominator of the fraction in 
the l.h.s. of Eq.iA.i9i 

Using Eqs. ([47|), ([43|), and ([48|) one can obtain that the 
energy of the compound solitori(pl3|), i.e. the soliton which 
belongs to the high-velocity branch, can be presented in 
the form 



vP> - Af>(u) 



(49) 



where 

P> = 
+ 



72 
1 



A-i 



/(7)-(2^-^7)7^(7) 



(50) 



is the canonical momentum of the compound soliton and 
1 



M>{v) = 



2880 



93320^3 



4wVl-27 (37^ -37 + 2) 
{A' -2A + A0){l-n{-f)) 



-A{A + 20)-/n{j) - 



iA-A)^A + 5) 



Hi) (51) 



is the corresponding nonlinear part. 

In the velocity interval v > Vc {A < 4) Eq. (^) for all 
values of the inverse radius a has a solution which tends to 
zero when d — > 00. There is an one-to-one correspondence 
between the velocity of soliton and its energy. But when 



(4c 



1) ^_6£ 
;2 - 1 - 729 



(52) 




1.60 



1.65 



1.70 



1.75 



Fig. 8. Energy of the pulse solitons versus velocity in the an- 
harmonic lattice with Kac-Baker (s — 0) long-range interaction 
calculated numerically (full line) and analytically (dashed line) 
for a = 0.05 and J = 0.0016. 



Eq. ( p8| ) has also two solutions in the interval v < Vc {A > 
4) . In other words under the condition (|5|) , which taking 
into account the explicit form of the speed of sound (^7|) , 
can be written as 



J > 



8{al ~ a2) 



(53) 



with ao ~ 0.16, there exist three stationary pulse solitons 
with different energy for each value of the velocity and the 
function H{v) has a characteristic Z-shape (see Fig. m. 
Comparison the dependence H{v) obtained analytically 
from Eqs. (|40|), (|4^)-(|5l|) with the dependence obtained 
numerically shows a good agreement between these two 
approaches (see Fig. |]). The function (|53| ) also describes 
very well the boundary which separates the area where 
H{v) has a iV-shape and a Z-shape (see Fig. 3 in Ref. 

0). 



4 Conclusion 

We investigated the effect of length scale competition in 
two systems which can bear pulse-like nonlinear excita- 
tions: a chain of coupled nonlinear oscillators (nonlinear 
Schrodinger model) and an anharmonic chain with cu- 
bic nearest-neighbor interaction. Two types of dispersive 
interactions were studied: the Kac-Baker interaction and 
screened Coulomb interaction. We have shown that the 
existence of two competing length scales: the radius of the 
dispersive interaction and lattice spacing provides the ex- 
istence of two types of pulse solitons one of which is broad, 
continuum- like and uni-component, and another is nar- 
row, discrete-like and has a compound structure. It con- 
sists of two components: short-ranged breather-like com- 
ponent which dominates near the center of the excitation 
and long-ranged component which contributes mainly to 
the tails of the excitation. These two stable solitons are 
separated by gap: an interval of spectral parameter (ve- 
locity in the case of Boussinesq equation and frequency in 
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the case of Nonlinear Schrodinger equation) where there 
are no stable stationary excitations. We show that the 
gap decreases when the radius of the dispersive interac- 
tion increases and it closes when the radius of the disper- 
sive interaction exceeds some critical value. In contrast 
to systems with short-range and moderate-range disper- 
sive interactions, where there is unambiguous correspon- 
dence between the velocity (frequency in the NLS equa- 
tion) of the soliton and its energy, in systems where the 
length scales differ significantly (ultra-long-range interac- 
tion) there is an interval where three different states with 
different energies (momentum, number of quanta, etc.) 
correspond to each value of the spectral parameter. This 
is a generic property of all systems with competing length 
scales and is due to the compound structure of excitations. 
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